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*This  report,  written  during  a  summer  appointment  at  MRC,  is  based  on  the 
Ph.D.  dissertation  of  C.  R.  Kennedy.  Tt  will  be  presented  at  a  conference 
on  dyneunical  systems  (Asilomar,  Dec.  9-14,  1979)  euid  published  as  a  joint 
paper  by  SIAM. 


SIGNIFICANCE  AND  EXPLANATION 


The  70 's  have  seen  a  surge  of  interest  in  the  behavior  of  solutions  of 
systems  of  differential  equations  which  in  some  sense  model  the  physics  and 
chemistry  of  diffusion  and  reaction.  One  line  of  work  has  grown  from  Turing's 
approach  to  morphogenesis,  which  he  saw  as  arising  from  stable  non-uniform 
solutions  in  an  otherwise  uniform  environment.  Another  was  stimulated  by  the 
Belousov- Zhabotinski  reaction,  a  colorful  oscillating  reaction  of  moderately 
exotic  chemicals,  and  by  the  discovery  that  catalytic  reactions  were  more  prone 
to  oscillation  than  had  been  realized.  Both  lines  are  obviously  connected  with 
bifurcation  theory  (the  latter  particularly  with  that  of  Poincare,  Andronov  and 
Hopf),  but  the  developing  understemding  of  strange  attractors  has  also  begun  to 
play  a  part. 

The  present  example,  basically  due  to  Rossler,  is  eui  attempt  to  exhibit 
this  behavior  in  as  simple  a  system  as  possible.  Two  equations  with  three  para¬ 
meters  represent  a  system  of  two  "reactions",  only  one  of  which  has  a  mild  non¬ 
linearity.  One  of  the  parameters,  corresponds  to  the  ratio  of  the  characteristic 
times  in  the  two  equations  and  has  no  effect  on  the  position  of  the  steady  state 
but  does  control  its  stability.  The  character  of  the  solution  in  one  cell, 
several  coupled  cells  or  a  continuum  under  variation  of  this  parameter  is 
studied;  it  is  found  that  the  uniform  steady  state  breaks  down  either  into  a 
uniform  oscillation  or  non-uniform  steady  state  amd  that  the  latter  further 
bifurcates  into  increasingly  complex  behavior. 


The  responsibility  for  the  wording  and  views  e 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


BIFURCATIONS  OF  A  MODEL  DIFFUSION- REACTION  SYSTEM 

R.  Aris* 


INTRODUCTION. 

The  bifurcation  of  solutions  to  equations  representing  the  interplay  of 
diffusion  and  reaction  has  received  a  great  deal  of  attention  in  recent  years 
both  from  a  theoretical  (1,  2,  3,  4)  and  from  an  e^qjerimental  (5,  6)  point  of 
view,  while,  since  the  discovery  of  chaos  (7,  8)  many  have  wondered  whether 
the  strange  flickerings  and  turbulent  behavior  of  quite  ordinary  systems  may 
not  be  ascribed  to  some  such  cause.  Rossler  has  demonstrated  how  several 
types  of  system  might  be  patient  of  strange  attractors  (9,  10)  and  we  wish 
to  pursue  one  of  his  suggestions  here.  Marek  has  considered  a  somewhat 
similar  system  (11)  while  a  bimolecular  reaction  has  engaged  Othmer  cind 
Kadas  (12). 

THE  SYSTEMS  CONSIDERED. 

The  vector  x  may  be  thought  of  as  a  vector  of  compositions  in 
with  components  x^,x^, . . .  ,x*'.  It  may  exist  in  three  contexts  which  are 
worth  comparing:  the  kinetic,  where  it  is  a  function  only  of  time,  x(t); 
the  compartmental ,  where  it  also  depends  on  an  enumerative  index,  Xj(t), 
j  =  1,...,N;  and  the  distributed,  where  (for  sinplicity)  it  depends  on  one 


position  variable,  z  say, 
guished  by  their  equations: 

as 

x(z,t).  Ihe  three  habitats  will  be 

distin- 

Kinetic  system: 

X  = 

f  (?) 

(K) 

Compartmental  system: 

X  .  = 

‘'D 

=  D(x  -  2x,  +  X  )  +  f(x.) 

-  ']+!  ~j  ~j-l  ~  ~j 

(C) 

-1'  ~N+1  "  ~N 

Distributed  system: 

X  = 

Dx  +  f(x),  X  =  0,  z  =  0,L  . 

~~zz  ~  ~  ~z 

(D) 

We  observe  that  all  three  systems  are  closed  except  in  so  far  as  the  function 
f  may  contain  terms  which  imply  the  removal  or  supply  of  any  species.  The 
kinetic  system  is  the  single  compeirtment;  the  distributed  is  the  limit  of  an 
infinite  number  of  compartments  if  D  is  made  proportional  to  (N/L) .  Of 
course  there  are  more  subtle  ways  of  approximating  the  distributed  system  by 
a  discrete  one  and  we  shall  use  a  Galerkin  method  in  calculations  on  D. 

The  equilibrium  or  steady  state  of  K  is  x^  given  by 

f(x  )  =  0  (1) 

'  ~e 

and  if  x  =  x  +  <fi  the  linearization  about  this  state  is 
-  -e  - 

=  Jifi  (2  ) 

where  J  is  the  Jar  bian  matrix  Vf(xg).  Local  stability  is  then  determined 
by  the  eigenvalues  ,f  J.  These  solutions  are  also  steady  states  of  (C)  and 
(D)  with  Xj  =  Xg,  J  =  i,...,N,  x(z,t)  =  Xg,  0  ^z  ^  L,  and  will  be  called 
uniform  steady  states.  If  tCj,  j  =  1,2,... ,N  are  eigenvalues  of  the  NxN 
matrix 
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then  those  of  the  linearization  of  the  uniform  state  of  C  are  given  by 

|Xl  +  -  J|  =  0,  j  =  1,2,. ..,N  (3) 

Since  the  <j  are  distinct  positive  numbers  with  =  0,  the  first  equation 
gives  the  eigenvalues  of  (K)  and  the  remaining  (N  -  l)n  eigenvalues  depend 
on  D.  Thus  the  stability  of  (K)  is  a  necessary,  but  not  sufficient,  condi¬ 
tion  for  the  stability  of  (C).  Moreover  it  is  differences  of  diffusivity  that 
destabilize  (C)  when  (K)  is  stable,  for  if  D  =  D1  the  nN  roots  of  (3)  are 

X,  “  tc.O,  1  ®  l,...,n,  3  “  1,...,N  , 

where  are  the  eigenvalues  of  J.  Similarly  the  linearization  about  the 

uniform  steady  state  of  (D)  is 

-e  IV  +  (4) 

~  -~zz 

and  a  perturbation  of  the  form  ce^^cos  (nirz/L)  is  stable  if 

I  XI  +  (mr/L)^D  -  j|  =  0  (5) 

has  roots  with  negative  real  parts  for  all  n.  It  follows  that  the  stability 
of  the  uniform  steady  state  of  (D)  bears  the  same  relation  to  (K)  as  does  (C). 

(C)  and  (D)  also  have  non-uniform  steady  states  unrelated  to  those  of  (K) . 
They  appear  in  symmetric  pairs,  for  if  satisfy 

0  =  D(x  -  2x  .  +  X  .  ,)  +  f(x  .)  (6) 

~  -e,3+l  'e,D  -e,3-l  -  -e,] 

then  so  do  x*  .  =  x  „  ,  and  if  x  (z)  satisfies 

'e,j  -e,N+l-j  -e 

0  =  D(x  )  +  f(x  )  (7) 

~  -e  zz  ~  ~e 

so  does  x*(z)  =  X  (L  -  z). 

~e  -e 

If  (K)  has  periodic  solutions  x(t)  =  g(t)  =  g (t  +  T) ,  then  (C)  and  (D) 
will  also  have  uniform  periodic  solutions  independent  of  j  and  z  respec¬ 
tively.  Similar  remarks  apply  to  the  linearized  stability  which  must  of 
course  be  determined  by  the  Ploquet  multipliers. 

ROSSLER'S  MODEL 

For  the  model  system  proposed  by  Rossler  (10)  we  will  take  x  =  (u,v) 
and  then  (K)  is 

u  =  u  +  1  -  uv/(k  +  u)  (8) 


V  =  v(ou  -  v) 


-2- 


and  has  three  parameters  k,v,  and  o.  On  the  first  and  last  of  these  the 
steady  state  depends  and  is  unique  and  "physical"  (i.e.  u  >  0,  v  >0) 
when  o  >  1 

V  =  au  =  {  (1  +  k)  +  /(I  +  k)^  +  4K(a  -  1) }  (10) 

e  e  2  (o  -  1) 

It  can  be  most  easily  found  as  the  intersection  of  the  line  v  =  au  with  the 
hyperbola  v  =  u  +  (1  +  k)  +  ku"-^  as  shown  in  Fig.  1. 


Of  the  conditions  for  stability  one  (that  the  product  of  the  eigenvalues 
be  positive)  is  automatically  satisfied  and  the  other  is  satisfied  when 

2 

V  >  g(K,a)  =  (u  -  k)/u  (<  +  u  )  (11) 

e  e  e 

where  Ug  is  given  by  (10).  The  dotted  lines  in  Fig.  1  are  contours  of  v 
and  in  Fig.  2  the  surface  g  is  shown.  If  the  parametric  point  lies  below 
the  surface  shown  in  Fig.  2  the  unique  steady  state  is  unstable.  Moreover  as 
we  cross  this  surface  the  eigenvalues  are 


+  ito  =  +  i/ (u^  -  k)  (2ic  +  u  +  ku  )/u  (u  +  tc)  (12) 

—  —  e  e  e  e  e 


Applying  the  criterion  for  stability  of  the  bifurcating  limit  cycle  shows  that 
it  is  stable  for  v  <  g(K,a)  and  in  the  neighborhood  of  that  surface.  Thus 
as  V  decreases  the  steady  state  loses  its  stability  euid  sheds  a  stable  limit 
cycle.  Fig.  3  shows  this  for  a  =  2.5,  <  =  0.144  for  which  Ug  =  0.87, 
ve  =  2.18  and  g  =  0.7;  limit  cycles  for  v  =  0.6,  0.5  and  0.3  are  shown. 
These  become  larger  and  more  triangular  as  v  ->•  0. 


STEADY  STATES  OF  ROSSLER' S  MODEL  IN  COMPARTMENTS. 


For  N  =  2  the  equations  are: 


Uj^  =  (0(0^  -  u^) 

+ 

1  .Ui 

-  +  Uj^) 

(13) 

Vi  =  ^(v^  -  V^) 

+ 

V(OUj^  - 

(14) 

U2  =  io(u^  -  u^) 

+ 

1.U2 

-  U2V2/«  +  u^) 

(15) 

^2  =  -  ^2’ 

+ 

v(au2  - 

"2> 

(16) 

For  the  moment  we  will  follow  Rossler  in  considering  only  the  exchange  of  the 
second  chemical  species  and  so  put  (2  =  0.  Since  (14)  and  (16)  are  linear 
they  can  be  solved  for  the  steady-state  v^^^  and  v^^  in  terms  of  u^^^  and 


V  . 

Ds 


a(tJ  t  v)  ^ 

V  +  2y  je 


+ 


oti 

V  +  2y  '^ie' 


i  ^  j;i»j  = 


(17) 


Substituting  into  (13)  and  (15)  then  gives 


For  the  uniform  steady  state,  the  characteristic  equation  has  two 
quadratic  factors: 

+  (v  -  g)A  +  v(ah  -  g)}{A^  +  (v  -  g  +  2p)  +  v(oh-g)-2)jg}  =  0.  (26) 

The  first  of  these  is  the  same  as  for  (K)  and  oh  -  g  >  0,  o  >  1.  Since 
p  >  0  we  now  have  the  conditions 

V > Max{g,2pg/ (oh  -  g) }  (27) 

for  stability. 

Again  we  expect  interesting  behavior  as  v  decreases  and  the  bifurcation 
picture  will  be  governed  by  the  order  in  which  the  coefficients  in  the  two 
quadratic  factors  change  sign.  There  are  two  possibilities  which  can  best  be 
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illustrated  by  Fig.  4  which  is  drawn  for  the  same  fixed  values  of  <  and  a 
as  Fig.  3.  The  exchange  parameter  y  is  also  held  constant  so  that  decreas¬ 
ing  V  is  represented  by  a  parametric  point  travelling  horizontally  to  the 
left.  The  uniform  steady  state  is  stable  to  the  right  of  ABC  and  if  the  para¬ 
metric  point  cross  AB  the  uniform  state  loses  its  stability  to  the  uniform 
periodic  state  in  which  the  two  compartments  behave  exactly  as  in  (K) .  The 
line  (0)BC  is  v  =  2yg/(oh-g)  so  that,  if  the  parametric  point  cross  BC,  the 
last  coefficient  in  (26)  changes  sign  first  and  the  system  loses  stability  by 
one  real  eigenvalue  becoming  positive  (rather  than  a  pair  of  complex  conjugates 
passing  across  the  imaginary  axis,  as  in  the  crossing  of  AB)  and  the  bifurcating 
solution  is  a  stable  non-uniform  steady  state.  In  the  region  DEBC  the  non- 
uniform  steady  states  are  stable  and  the  same  state  obtains  over  any  ray  through 
the  origin,  since  u+  and  u_  depend  only  on  v/y,  as  shown  in  Fig.  5.  They 
lose  their  stability  across  the  curves  DE  and  EB  as  will  be  discussed  after  a 
look  at  the  steady  states  of  the  system  with  N  >  2 . 


For  N=3,  a  quadratic  equation  can  be  found  for  U2e  (^t  the  ejqiense 
of  considerable  algebra)  and  further  equations  for  U]^g,  113^.  They  give  rise 
to  three  classes  of  steady  solution: 

I.  uniform  steady  state  (u  ,u  ,u  ) 

e  e  e 

II.  asymmetric  steady  states  (u,  ,u„  ,u,  ),  (u,  ,u.,  ,u,  )  ,  u,  ^  u, 

le  2e  3e  3e  2e  le  le  3e 

III.  symmetric  steady  states  (u,  ,u^  ,u,  ) ,  u,  =  u, 

le  2e  3e  le  3e 

Vte  shall  not  give  the  formulae  but  show  in  Fig.  6  the  solutions  that  make 
all  the  Ujg  positive.  For  this  set  of  parameters  only  Class  I  can  obtain 
for  v/y  >  1.45;  I  eind  III  are  allowed  for  .57  <  v/y  <  1.45  and  all  three 
are  possible  for  v  <  .57y. 

For  large  N  we  can  simplify  the  picture  by  being  more  general  and 
allowing  the  exchange  of  both  components  i.e.  C  ^  0.  Then,  by  (3),  the 
stability  of  the  uniform  steady  state  is  governed  lay  N  quadratics,  j  =  l,...,N, 

2 

X  +  {v  -  g  +  K .  (u  +  y)  }X  +  v(ah  -  g)  +  k.  (ioyic.  +  uv  -  yg)  =  0 


Since  =  0  and  the  k j ,  j  >1,  are  positive  the  coefficient  of  X  changes 
sign  first  for  j  =  1  as  v  becomes  less  than  g.  If  all  the  constant  terms 
are  positive,  as  will  be  the  case  when  S)  eind  y  are  sufficiently  small,  the 
imiform  steady  state  gives  rise  to  a  uniform  periodic  one.  Across  the  line 

V-  -  2 

ah  -  q  +  u)K . 

—  _ 1 

K .  (g  -  (jjK . ) 

3  3 

there  is  a  bifurcation  to  a  non-uniform  steady  state.  (Min^  is  the  least 
positive  value.)  Since  we  are  primarily  interested  in  the  dynamic  behavior  we 
will  not  try  to  explore  this  complex  picture,  but  choose  oi  so  as  to  simplify 
it.  We  note  that  j  =  2  may  not  be  the  minimizing  value  of  j  in  (28).  For 
example  with  N  =  10  amd  u  =  2  only  one  of  the  expressions  in  (28)  is  posi¬ 
tive  and  we  have  a  diagram  like  that  in  Fig.  4.  The  non-uniform  steady  states 
in  the  DEBC  region  are  the  same  along  rays  of  constant  v/y.  Three  exaitples 
are  given  in  Fig.  7. 


V  Min  M. 
3>1  ' 


(28) 
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Figure  4  Bifurcation  diagram  in  the  vi,v-plane  for  M=2,  ic- 0.144,  o 


03  0'95 


I  Figure  5  Uni'"<.  and  non-uniform  steady  states  as  functions  of  v/p. 

I 


I 

i 

I 


Uniform  and  non-uniform  steady  states  for  N  =  3.  I  denotes  the 
uniform  steady  state,  II  denotes  the  axisymmetric  steady  state,  III 
the  symmetric  one. 


Figure 


CONCENTATION  V:  CONCENTRATION 


DYNAMICS  OF  ROSSLER  IN  COMPARTMENTS. 

Rossler  gave  several  examples  of  the  dynamics  with  N  =  2  and  we  shall 
follow  his  lead  by  systematically  varying  v  for  N  =  2  before  increasing 
N.  In  particular  we  keep  k  =  0.144,  o  =  p  =  2.5  and  u  =  0,  so  that  we 

move  to  the  left  on  the  line  (L)  of  Fig.  4.  To  the  right  of  F  there  is  only 

the  uniform  state  U]^  =  U2  =  0.87,  Vj  =  V2  =  2.18.  Between  F  and  G  the 

non-uniform  steady  states  can  be  followed  in  Fig.  5,  but  at  g  this  becomes 
unstable  by  a  Hopf  bifurcation,  p  being  0.755.  The  kind  of  oscillation 
that  arises  is  shown  in  Pig.  8  for  which  v  =  .75  (point  H  in  Fig.  4).  At 
V  =  .696  (point  J  of  Fig.  4),  this  oscillation  loses  its  stability  as  a 
Floquet  multiplier  passes  out  of  the  unit  circle  through  -1,  and  a  stable 
periodic  solution  of  twice  the  period  arises;  this  solution,  for  v  =  0.690, 
is  shown  in  Fig.  9.  By  tracking  the  critical  Floquet  multiplier  as  in  Fig.  10, 
we  find  a  bifurcation  which  again  doubles  the  period  to  about  four  times  the 
original  one.  The  stability  of  this  quickly  breaks  down  and  the  period  again 
doubles. 

At  this  point  it  becomes  difficult  to  keep  track  of  what  is  happening  for, 

even  when  there  is  an  attracting  orbit,  it  takes  an  unconscionably  long  time 

to  lock  onto  it.  Vfe  tried  recording  the  values  of  everything  each  time  u^^ 
passed  through  a  maximum  eind  this  poor-man' s- Poincare -map  gives  clear  support 
for  the  picture  just  drawn.  There  was  less  clear  evidence  of  period  3  oscil¬ 
lations  but  much  indication  of  irregular  behavior,  though  whether  this  is 
chaos  or  a  solution  of  long  period  cannot  be  determined  empirically.  It  may 
be  that  as  v  decreases  we  have  a  sequence  of  bifurcations  such  as  is  known 
(12)  to  obtain  for  Xn+1  =  Xxn  d  -  Xn)  or  there  may  be  a  Hopf  bifurcation  of 
the  Poincare  map  to  a  torus  eind  so  on.  It  is  also  possible  that  the  successive 

bifurcations  cease  altogether  for,  when  v  <  0.515  (point  K  in  Fig.  4),  the 

uniform  periodic  solutions  are  stable.  There  would  thus  appear  to  be  a  locus 
EKM  in  Fig.  4  such  that  the  uniform  periodic  solution  is  stable  in  the  region 
OABEKM.  The  structure  of  left  hand  part  of  the  zone  MKEGD  is  as  complex  as 
it  is  mathematically  erogenous. 

Explorations  for  N  =  10  give  evidence  that  the  same  general  behavior  is 
to  be  expected.  In  the  case  k  =  0.144,  vi  =  250,  u)  =  2 ,  o  =  2.5  the  stable 
non-viniform  limit  cycles  in  five  cells  of  the  symmetric  system  are  shown 
in  Fig.  11  for  v  =  0.6.  By  v  =  0.57  the  solution  has  already  bifurcated 
twice.  For  yet  smaller  values  of  v  the  solution  is  shown  in  Fig.  12,  where 
the  smaller  ripples  (only  uj^  is  shown)  are  'near'  the  uniform  periodic  solution 
and  the  large  departures  'near'  the  non-uniform.  Since  the  uniform  periodic 
solution  is  ultimately  stable  the  scheme  of  things  would  seem  to  be  as  sketched  in 
Fig.  13,  where  as  v  decreases  the  increasing  stability  and  attractiveness  of 
the  uniform  orbits  tears  the  non-uniform  orbit  apart.  This  is  in  keeping  with 
Rossler' s  observations  (10). 


ROSSLER  DISTRIBUTED. 

Here  we  will  take  C  =  z/L  and  work  with 

u  =il)u_.  +  l+  u-  uv/  ( ic  +  u ) 
I 

V  =  yv  +  v(cru  -  v) 

T 


(29) 

(30) 
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Figure  9  Non-uniform  oscillations  of  period  -2T  in  the  two  cells. 


Figure  10  Critical  Flcxjuet  multiplier  and  period  as  functions  of 
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unifom  pericxiic  solutions  in  five  of  the  ten  cells  shown  as  phase-plane  diagrams 


First  Cell  Concentration  Ui  First  Cell  Concentration 


T- Dimensionless  Time 

Figure  12  Time  variation  of  the  concentration  uj^  for  em  orbit  which  spends 
most  of  its  time  near  the  uniform  periodic  solution. 
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Figure  i:?  Schematic  of  how  the  attracting  orbits  evolve. 


with 


u^=v^=0  at  C=0,1 


(31) 


The  uniform  solution  Ue(C)  =  Ue»  Vg(C)  =  Vg  is  still  given  by  (10),  while  by 
(5),  its  stability  is  governed  by 

2  2  2  2  2  2  2 

A  +  (v  -  g  +  ((jj  +  v)n  TT  )  A  +  [v  (ah  -  g)  +  n  ir  (n  it  uy  +  £ov  -  yg)  ]  =  0  (32  ) 


If  (j  and  y  are  non-zero,  these  eigenvalues  tend  to  -n  it  Si  and  -n  it  y 
(i.e.  to  -“>)  as  n  but  if  fij  =  0  one  tends  to  g  and  the  other  to 

+00.  We  shall  avoid  the  case  of  £j  =  0  for  this  reason. 

We  again  as)c  how  the  uniform  solution  bifurcates  and  again  have  a  diagram 
as  in  Fig.  4  for  which  the  line  OBC  is 


Min’*” 

^  2  2-1 
ah  -  g  +  n  n  (J 

Min'*'  M(n) 

/  ~  2  2,  2  2 

V  — 

n>0 

(g  -  (ijn  TT  )n  TT 

n>0 

(32) 


Since  g  <  1,  a  necessary  condition  for  instability  by  violation  of  (32)  is 
uiir^  <1.  If  y  is  small  enough,  reduction  of  v  will  take  the  solutions  of 
the  quadratic  for  n  =  0  into  the  unstzdale  half  plane  as  a  pair  of  complex 
conjugates.  If  y  is  above  the  point  B,  a  disturbance  of  wave  number  (nQ 
say)  that  minimizes  the  right  hand  side  of  (31)  begins  to  grown.  This  will  be 
illustrated  by  some  numerical  examples,  but  we  mention  in  passing  that,  using 
the  usual  type  of  analvsis,  a  criterion  to  distinguish  between  sub-  and  super¬ 
criticality  was  obtained.  It  is  too  complicated  a  formula  to  permit  a  general 
conclusion,  but  it  is  suspected  that  as  v  decreases  the  uniform  steady-state 
gives  up  its  stability  to  the  non-uniform  and  that,  whilst  the  uniform  state  is  stable, 
stable  non-uni  form  states  do  not  exist  when  Uq  is  odd  but  do  exist  when  nQ  is  even . 

Three  cases  were  calculated,  with  a  =  2.5,  ic  =  0.144  amd  Si  =  0.01, 

0.02  eind  0.05.  When  (5=0. 01,  0<M(2)<M(1),  M(n)<0,  n>2,  so  that 
the  minimizing  ng  is  2  and  we  should  e3q)ect  a  symmetric  solution  to  arise 
first  and  perhaps  give  its  stability  to  am  unsyrametric  one  as  v  decreases 
further  since  M(l)  is  also  positive.  Fig.  14  is  the  result  of  transient 
calculations  which  converged  to  the  stadale  non-uniform  steady  states.  Though 
these  are  far  from  the  transition  it  is  the  larger  v  that  gives  the  sym¬ 
metric  state  (n  =  2).  There  is  an  analogy  with  the  case  of  three  compart¬ 
ments  where,  as  v  decreases,  the  symmetric  solutions  (class  III)  are 
superceeded  by  the  unsymmetric  (class  II). 

For  0)  =  0.02,  ng  =  1  amd  M(n)  <0  for  n  >  1.  The  non-uniform 
steady  states  depend  only  on  the  ratio  y/v  and  three  are  shown  in  Fig.  15 
that  correspond  to  the  rays  in  Fig.  16.  The  eigenvalues  of  the  linearization 
about  the  non-uniform  steady-state  were  calculated  and  it  was  found  that 
stability  was  lost,  by  the  real  part  of  a  pair  of  complex  conjugate  eigen¬ 
values  becoming  positive,  on  the  line  BED.  Since  it  was  confirmed  that  the 
uniform  periodic  solution  was  stable  for  small  enough  v,  it  seems  likely 
that  there  is  a  region  of  complex  behavior,  KED,  just  as  in  the  case  N  =  2. 
Calculations  with  Si  =  0.05  concur  2md  the  agreement  between  calculations 
for  N  =  10  euid  (D)  give  some  confidence  that  this  case  of  (C)  gives  a  good 
indication  of  the  dynamics  of  (D) .  Ceitainly  the  worst  can  be  expected  of 
(D)  for  which  Pismen  (13)  has  shown  how  chaos  may  arise. 
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Symmetric  and  asymmetric  steady-states  for  the  distributed  system 


Figure  16  The  presumed  bifurcation  diagram  for  the  distributed  system.  T 
rays  labeled  i,  2,  and  3  are  the  loci  of  the  solutions  shown  in 
Figure  15.  The  uniform  solution  is  stable  to  the  right  of  ABC, 
the  uniform  periodic  solution  to  the  left  of  KABA.  The  non- 
uniform  steady  solutions  are  stable  within  the  arms  of  DEBC  and 
a  region  of  exotic  behavior  exists  in  KED. 
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